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ABSTRACT 

The interpretation of gravitational microlensing events caused by planetary 
systems or multiple stars is based on the n-point-mass lens model. The first plan¬ 
ets detected by microlensing were well described by the two-point-mass model of 
a star with one planet. By the end of 2014, four events involving three-point- 
mass lenses had been announced. Two of the lenses were stars with two planetary 
companions each; two were binary stars with a planet orbiting one component. 
While the two-point-mass model is well understood, the same cannot be said for 
lenses with three or more components. Even the range of possible critical-curve 
topologies and caustic geometries of the three-point-mass lens remains unknown. 

In this paper we provide new tools for mapping the critical-curve topology and 
caustic cusp number in the parameter space of n-point-mass lenses. We perform 
our analysis in the image plane of the lens. We show that all contours of the 
Jacobian are critical curves of re-scaled versions of the lens configuration. Uti¬ 
lizing this property further, we introduce the cusp curve to identify cusp-image 
positions on all contours simultaneously. In order to track cusp-number changes 
in caustic metamorphoses, we dehne the morph curve, which pinpoints the posi¬ 
tions of metamorphosis-point images along the cusp curve. We demonstrate the 
usage of both curves on simple two- and three-point-mass lens examples. For 
the three simplest caustic metamorphoses we illustrate the local structure of the 
image and source planes. 

Subject headings: gravitational lensing: micro — planetary systems — stars: 
multiple 
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Introduction 


The effect of the gravitational field of an astrophysical object on the passage o f light 


from a backgrou nd source was first computed for a point-mass-like object flEinsteinl Il936 


Tjnklll936l. 119371) . which provided the simplest model of a gravitational lens. The model can 
be used in some cases as a first approximation for describing the strong lensing of a quasar 
by a galaxy. Its main application, however, has been in Galactic gravitational microlensing, 
where it accurat ely describes the lensing of a background star by another star passing in the 
foreground fe.g.. IPaczvhskil 1199611 . 


The more than ten thousand observed microlensing events 

nclude several hundred 

cases of lensing by binary stars (lists can be founc 

in i 

Vlcock et ah 

2000; 

Jaroszvhski 

2002; 

Jaroszvhski et ah 

2004, 

2006, 

2010; 

Skowron et ah 

200 

7). In addition, there are more than 


30 cases of microlensing by a star with a planet (e.g.. lBond et al.ll2004 iBeaulieu et al.ll2006l) . 


Planets are currently the primary goal of microlensing surveys, which are sensitive to Solar- 
system analogues including Earth-mass planets in AU-scale orbits. Binary stars and stars 
with a planet can be both described by two-point-mass microlensing, which has been well 
understood theoretically for the past twenty years. 


Chang fc Refsdall (1l984j) first described an extension of the single point-mass lens ob¬ 
tained by adding a constant external shear. Such a syste m has extended caustics with 


different geometries in the weak and strong shear regimes. ISchn eider fc WeissI fll986l) per¬ 
formed a detailed analysis of a two-point-mass lens with equal masses. They demonstrated 
that there were three distinct regimes for the critical curve and caustic: close, intermediate, 
and wide (in order of increasing separation of binary components). Tr ansitions between 


neighb oring regimes involved beak-to-beak metamorphoses of the caustic. lErdl fc Schneider 


(I 1993 I) studied the lensing by two independent objects generally in different lens planes. As a 
special case they included the analysis of the single-plane (thin) two-point-mass lens with an 
arbitrary mass ratio, showing that the regimes ar id metamorphos e s wer e the same as in the 
equal-mass case. This result was corroborated bv lWitt fc PettersI (1199311 . From the range of 
subsequent papers we mention here the analysis by DominilJ ~f 1999 ) of special cases of the 
binary lens. 

However, alr eady at n = 3 obse rvations overtoo k theory. Three-c omponent lensing was 
detected first by Gaudi et ah ( 2008 ) and second by Han et aP J2013). B oth cases involved 


micro lensing by a star with two planets. Most recently, iGould et al.l (120141) and iPoleski et ah 


(120141) documented two events with microlensing by a binary star with a planet orbiting one 
of the components. Despite these observational advances, we still have no general insight into 
the possible regimes, types of critical curves and caustics permitted by the three-point-mass 
lens model. Nevertheless, a number of theoretical papers have been published on different 
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aspects of triple-lens microlensing. 


The first step beyond two-point-mass leasing was taken bv iGrieger et al.l fjl989[) . who 
studied the effect of an additional constant shear and demonstrated the appearance of the 
swallo w-tail and butterfly ca ustic metamorphoses. This particular model was studied late r 
also bv lWitt fc PettersI (119931) . The general triple-lens equation was discussed bv iRhiel (1200211 . 
Other works concentrated on specihc types of triple lens es. Most frequent among them are 
papers studying lensing by a star with two p anets (e.g., iGaudi et al.lll998l: iHan et al.l 12001 


Rvu et ahl I 2 OIII: ISong et al.l |2014J: IZhu et ah 
binary star with a single planet (e.g.. iLee et al 


20141). Several pamp ers explored lensing by a 


20081: lHanll2008al) . Finally, the possibility of 


detecting ex o-moons in microle n sing events in volvi ng a star with a planet with a moon was 
discussed by Han &: Han ( 2002I) . Han (2008b), and Liebig fc Wambsganss ( 2010 ). 


Much of the literature on general n-point-mass lenses revolved around the question of 


the maximu m poss i ble number of images ge n eratec 


jy su c h a 


ens. Individual contribu- 


tions include 

R] 

Qie (199' 

7). Mao et al. 

(1997). R] 

tiie (2001. 

2003) 

, Khavinson & Neumann 

(200 

1). Asada 

l( 

2009bl). 

Luce et al. (S 

i014), and 

Sete et al. 

(201. 

5). Pivotal among these 

are 

Hhiel ( 2 OO 3 I). who coniectured th 

e maximum number of images was 5 (n — 1), and 
'oved the conjecture. 

Khavinson & Neumann 

2000). who p] 


Other papers u sed perturbat ive methods to study the n-point-mass lens f ormed by a star 
with n — 1 planets. iBozzal (119991 ) explored the caustics of such a system, and lAsadal (j2009al ) 
pres ented perturba tive solutions of the lens equation. Finally, for the general n-point-mass 
lens iBozzal (j2000al Jbl) studied the geometry of caustics in the close and wide limits. 


The initial goal of the research behind the present paper was to syst ematically explore 
the pr opert ies of three-poin t -mass lenses. Following the methods used by lErdl fc Schneider 
(II 993 I) and IWitt fc PettersI (119931 ) for the general two-point-mass lens, we set out to map 
critical-curve topologies and caustic geome tries in the p a rameter space of the t hree-p oint- 
mass lens. Some ear l y resu lts appeared in Danek ( 2010 ). Danek fc Hevrovskv (j2011 ). and 


Danek fc Hevrovskvl (120141) . In the course of our work we devised an efficient method for 
tracking the changes in critical-curve topology and caustic cusp number in the image plane 
of an arbitrary n-point-mass lens. We present these results here, while their application to 
the analysis of simple three -point-mass lens models appears separately in a companion paper 
( Danek fc Hevrovskv 2015 ). 


We start in § [2] by describing the studied lens model, a thin n-point-mass gravitational 
lens with no external shear and no convergence due to con tinuous matter. We transform the 
lens equation to complex formalism following (lWittlll990f) . which translates the analysis of 
lens properties to the study of polynomials, curves, and real and complex functions in the 
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complex plane. 

Our approach is based on investigating the properties of the Jacobian of the lens, viewed 
in §[3] as a surface (real function) over the complex image plane. The main new insight is the 
correspondence of Jacobian contours with critical curves of zoomed-in or zoomed-out lens 
configurations, as described in § 13.31 The nature of the critical curve, and the changes of its 
topology in the parameter space of the lens are discussed in § HI 

In §|3]we study the caustic of the lens, paying particular attention to the changes in cusp 
number during caustic metamorphoses. We introduce two further tools for finding important 
points along Jacobian contours in the image plane: the cusp curve in 15.31 for tracking the 
positions of cusp images, and the morph curve in 11 15.41 for pinpointing the images of caustic- 
metamorphosis points. In § 15.51 we give specific conditions for distinguishing butterfly and 
swallow-tail metamorphoses. 

We demonstrate the introduced method in § |6] with the example of a binary lens, and 
two triple lenses with swallow-tail and butterfly metamorphoses, respectively. Finally, we 
illustrate the local structures in the image and source planes in the vicinity of the three 
simplest caustic metamorphoses of the n-point-mass lens. The main results are summarized 
in § El 


2. Lensing by an n-point-mass lens 

In this work we study a composite gravitational lens formed by an association of n 
objects, which we approximate by point masses with masses Mj, j = 1.. .n. We retain the 
usual thin-lens approximation, in which variations in the line-of-sight distance of individual 
objects from the observer are neglected. The entire lens is thus assumed to lie at a single 
distance Di from the observer and Dig from the background source, so that Dg = Di + Dig 
is the distance from the observer to the source. The angular positions Oj of the individual 
objects (lens components) in the plane of the sky, the position 9 of an image formed by the 
lens, and the position (3 of the source are all defined relative to a specific point of the lens. 
The usual choices for this point are: the position of a lens component, the center of mass 
of the lens, or its geometric center. While the equations in this work are independent of 
the choice, in the presented examples we use the geometric center as the origin. All object 
positions must be distinct, 0j ^ Oj. for all j ^ k, in order to qualify as an n-point-mass lens. 
In the thin-lens approximation any pair of objects exactly aligned with the observer acts as 
a single object with the total mass of the pair, reducing the number of lens components by 


one. 
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The lens defines a characteristic scale in the plane of the sky: the angular Einstein 
radius corresponding to the total mass of the lens. It is given by 


Or — 


AGMtot Dis 


Di 


( 1 ) 


where G is the gravitational constant, c is the speed of light, and Mtot = YTj=i the total 
mass of the lens. In the absence of external shear and convergence due to continuous matter, 
the relation between the source and image positions is expressed by the lens equation 


^ a _ n , 

/3 = hi ^ 

i=i 


e-9,\^ 


( 2 ) 


Here fij = Mj/Mtot are the fractional masses in units of total mass (hence X]i=i hi = !)• 

In order to stud y the prope rties of the lens it is advantageous to switch from vectorial 
to complex notation fjWittI 1199011 . in which all angular positions in the plane of the sky are 
expressed in terms of complex variables. If we additionally measure all angular positions 
in units of Or, the positions of the individual components are zj = {Oji + i9jii)/9E, the 
image position 2 ; = {9j + i9ii)/9E, and the source position ( = {(3i + i(3 ii)/9e, where we 
used Roman indices to mark vector components, and i is the imaginary unit. The inverse 
transformation from complex to vectorial notation is obtained by combining a position with 
its complex conjugate. For example, the image position 6 = {z + z, —iz + iz) 9El‘i-i where 
2 denotes the complex conjugate of 2 ;. In matrix notation the transformation 6 ^ z can be 
described by (z, z) = G {Oj, 9n) with 


G = 


9, 


1 i 
1 -i 


( 3 ) 


The same matrix transforms similarly Oj Zj and j3 ^ The opposite transformation 
z ^ 6 can be described using the inverse matrix {9j, On) = G~^ (z, z), given explicitly by 


c-i = 


9f 


1 
— i 


Matrix G ^ transforms also Zj —)■ Oj and C —?■ /3 in an analogous wajj^. This descripti 
be useful for studying the Jacobian of the lens mapping and its properties in § [3 

The lens equation (|2]) attains a simpler form in complex formalism. 


( 4 ) 


ion will 






j=i j 


( 5 ) 


^In case the positions are already in units of 9e^ set 0^; = 1 in equations m and 
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The positions of images for a given sonrce position are usnally fonnd by transforming equa¬ 
tion ([5]) to polynomial form purely in terms of This can be achieved by taking the complex 
conjugate of equation ([S]), expressing ^ in terms of 2 ; and ( from it, and substituting back 
into equation ([3]), obtaining 


< = ^-T.7Z 


fXj 


;=1 ^ l^k=l z- 


l^k 

Z-Zk 


( 6 ) 


The expression on the right-hand side can be converted to a rational function of ^ by gradual 
transformation to a common denominator. Multiplying by the hnal denominat or leads to a 
polynomial equation of degree -f 1 for the image position z (e.g., IWittlll990h . 

Such an equation always yields a set of -|- 1 solutions. However, not all of these are 
true images. Additional solutions were introduced by the substitution for 0 when deriving 
equation ([S]), which is not an equivalent step. Only those solutions that satisfy the original 
lens equation ([5]) correspond to images of the lensed sour ce. The actual number of images for 
n > 1 ranges from n-l-ltohn — Sin increments of two flRhiell2003l: iKhavinson &: Neumann 
20061 ) . 


The number of images depends on the position of the source with respect to the caustic , 
which is a characteristic curve of the lens in the source plane (e.g., ISchneider fc Weisslll986l) . 
When a point source enters the caustic, a pair of images appears and branches off from a 
point on the critical curve of the lens in the image plane. Conversely, when a point source 
exits the caustic, a pair of images meets and vanishes at the critical curve. In special cases, 
when the point source enters or exits the caustic at an intersection point of two or more 
sections of the caustic, two or more pairs of images may appear or vanish simultaneously at 
different points along the critical curve. 


From the description above, it follows that the critical curve Zcc is a degenerate image 
of the caustic Cc- This property is usually used in reverse to hnd the caustic of the lens by 
setting z = Zcc in equation (|5]). In mathematical terms, the critical curve is defined as the 
zero set of the Jacobian determinant corresponding to the lens equation ([5]). 

For a given image of a point source its flux amplification is given by the reciprocal abso¬ 
lute value of the Jacobian determinant, and its parity is given by the sign of the determinant. 
Since the Jacobian determinant plays a pivotal role in the analysis of gravitational lensing, 
we explore its properties in more detail in the following section. 
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3. The Jacobian and its properties 
3.1. Definition and global character 


The Jacobi matrix in real vectorial leasing notation is defined as 

j ^ d {(3 1 , /3jj) 

d{ei,eu)' 

The usual definition of the Jacobi matrix in complex notation is 

- d{zrz) ’ 


(7) 


( 8 ) 


where z and z are taken as independent variables flWittlll990f) . The two matrices are con¬ 
nected by the simple transformation 


J = CJaC 


-1 


(9) 


with the matrices C and C~^ defined in equations ([3]) and (0]), respectively. Since equation (l9i) 
conserves the determinant, det J = det Jo , all properties of det Jo{0) as a real function over 
the real plane are identical to those of det J {z) as a real function over the complex plane. 
The determinant (hereafter the “Jacobian”) can thus be viewed and explored as a simple 
two-dimensional surface. 


In the case of our studied lens equation ([5]), the Jacobian simplifies to 


det J (z) = diQ - SjC 9--C = 1 - |9--Cr = 1 - 


i: 

i=i 


fXj 


( 10 ) 


where we used abbre viated notation for partial derivatives, e.g., dzC = d(/dz flWittlll990 


Witt fc Petterslll993l) . From the final expression we immediately see that det J {z) < 1 over 
the entire complex plane. The Jacobian is a continuous smooth function everywhere except 
at n poles at the positions of the lens components, where det J {zj) = —cxo. Far from all 
the lenses, for \z — Zj\ S> 1, it approaches unity, det J {z) ^ 1. The set of points with zero 
Jacobian forms the critical curve of the lens, det J (zcc) = 0. 


3.2. Stationary points: maxima and saddles 

The stationary points of the Jacobian fulfill the condition dz det J = 0. Since the 
Jacobian is a real function, dz det J = 0 is then implied automatically. Using equations ffTOj) 
















and ([5]) we obtain the expression 


dz det J = —dzC dzC = 0 , 


( 11 ) 


where we nsed abbreviated notation for the second derivative, = d‘^(/dz‘^. Eqnation ffTTl) 
yields two types of solntions, d^C = 0 and d‘^( = 0. The natnre of these solutions can be 
found from the Hessian matrix of the second derivatives along the real and imaginary axes, 
with components hij = d'^det Jo/dOidOj where i,j = I, II. Expressed in terms of complex 
derivatives, the Hessian matrix 


f d\ det J dl^deiJ \ 


( 12 ) 


where is the transpose of matrix C. Substituting for the Jacobian from equation ffTOj) . 
we obtain the Hessian determinant 

deti/ = T(|SJ^|‘_|4<;|2|Sj^p) , (13) 

E 

Note that the Hessian has the opposite sign when compared with the determinant of the 
central matrix on the right-hand side of equation flT^ . 

The stationary points with d^C = 0 have detiJ = > 0, hence they are 

local extrema if 7^ 0- From equation ffTOj) we see that all points with dzC = 0 have 
det J = 1. Therefore, these points are always maxima of the Jacobian irrespective of the 
second derivative, achieving the same peak value as for z —)■ cxo. Note that the Jacobian 
value also implies that these maxima never he on the critical curve. In terms of the lens 
parameters, the condition for Jacobian maxima is 


ax" = E r-^ 


/ij 




= 0 , 


(14) 


where we used the complex conjugate of equation ([S]). Multiplication by all denominators 
converts equation fllTll to a polynomial of degree 2 n — 2 with highest-degree-term coefficient 
equal to 1. The n-point-mass lens Jacobian thus has at most 2n — 2 different local maxima, 
since the polynomial generally may have multiple roots, for which = 0. In the 

limit when all lens components are mutually much closer than the Einstein radius, \zj\ 1, 
the Jacobian maxima lie at the centers of the small loops of the critical curve. 


The stati onary points with dIC = 0 have detiJ = —d6l'^\dz C\‘^\d'^C\‘^ < 0, hence they are 
saddle points (jSchneider fc Weisslll986l: lErdl fc Schneideiill993l) unless dzC = 0 or d^C, = 0. 
In the former case they are higher-order maxima, as shown above. In the latter case they 
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are higher-order saddle points, such as the monkey saddle if d^C, 7 ^ 0. In terms of the lens 
parameters the condition for Jacobian saddle points can be written as 



«= 5 : 

i=i 



0 , 


(15) 


where we used equation ([5|) and complex conjugation. Multiplication by all denominators 
converts equation flT^ to a polynomial of degree 3 n — 3 with highest-degree-term coefficient 
equal to 1. The n-point-mass lens Jacobian thus has at most 3n — 3 different saddle points. 
Some solutions of equation ffTSj) may correspond to higher-order maxima, or several solutions 
may correspond to a higher-order saddle. The saddle points play a key role in determining 
changes in the topology of the critical curve, as discussed further in § HI 


3.3. Jacobian contours and critical curves 


The n-point-mass-lens Jacobian has a remarkable property we have not found mentioned 
in the literature on gravitational lensing. We demonstrate here that all level curves of the 
Jacobian, seen as contours in the image plane, correspond to critical curves of re-scaled 
lens configurations. We can show this by taking any level curve of the Jacobian surface, e.g., 
det J (z) = A, where the Jacobian range restricts the parameter A to real values A < 1. Using 
the last expression for the Jacobian from equation fllOl) . the equation for the corresponding 
contour in the image plane is 


1 - 


n 



A. 


(16) 


If we subtract A and divide by 1 — A , we get 


1 - 


1 

1-A 


n 



0 . 


(17) 


This step can be taken anywhere except exactly at maxima. There A = 1 and the cor¬ 
responding contours are point-like. Now if we re-scale the image plane including the lens 
positions by a real positive factor. 


z' = z — A , Zj = Zj — A , 


equation flTTl) simplihes to 


1 - 


7Z1 (*' - 


(18) 


2 


0 , 


(19) 
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which is the critical-curve equation for a lens with components at z'y 

One can proceed similarly the other way around: if we change the scale of a lens 
conhguration by shifting all lens components from Zj to azj, where a is a real positive 
factor, we find that the critical curve is equal to the Jacobian contour det J = 1 — of the 
original configuration re-scaled by a. 

This property is especially useful for studies of the lensing behavior, and it further 
enhances the significance of the Jacobian surface. For any particular n-point-mass lens the 
surface includes as its contours not only the critical curve of the given conhguration, but 
also the critical curves of all directly similar conhgurations, arbitrarily shrunk or expanded. 
For values A —)■ — cxd in equation flT6|) we get the critical curve of a set of n individual widely 
separated lenses (the “wide” limit), while for A —)■ 1“ we get the critical curve of a combined 
lens with all components approaching a single point (the “close” limit), as can be seen from 
equation (lT8l) . 

In Figured] we demonstrate the Jacobian contour o critical curve correspondence on 
the example of a two-point-mass lens (hereafter for simplicity “binary” lens) with masses 
/xi = 4/5 and ^2 = 1/5 separated by sq = 1 Einstein radius. In the left panel we plotted 
the Jacobian surface det J( 2 ;), with its “trunks” around the lens positions truncated at 
det J = —15. On the surface we marked a sequence of Jacobian contours specihed in the 
caption. The highest contours outside the upper self-intersecting det J = 0.711 contour each 
have two additional small loops around the positions of the two Jacobian maxima. The three 
saddle points of the Jacobian surface lie at the self-intersection points of this contour and of 
the lowermost det J = —11.0 contour. 

In the right panel we plotted the same contours as seen from above in the image plane, 
adding the innermost det J = —100 curve that would he below the plot in the left panel. 
For illustration, selected curves corresponding to Jacobian values marked in the caption in 
boldface are labeled by the re-scaled binary separation sq ^1 — det J. The given contour is 
identically equal to the critical curve in a similarly re-scaled image plane. 

The right panel may also be viewed as a sequence of critical-curve plots for lenses 
with the same mass ratio and different separations s, with the position in the image plane 
measured in units of s instead of 6*^. Finally, the sequence of curves also corresponds to 
holding the absolute angular separation s9e fixed and changing the angular Einstein radius 
9e (e.g., by varying the total mass of the lens Mtot-, or any of the distances Dg, Di, Dig). 

Note that the positions of Jacobian maxima and saddle points scale together with the 
lens positions, i.e., their relative conhguration stays hxed with the lens components. This 
can be seen from equations ffTT|) and flT^ . in which any scaling factor in the denominator 
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drops out due to the zero right-hand side, so that both equations retain their form. The 
correspondence with contours of the Jacobian surface in the left panel of Figure [1] provides 
the best illustration: the maxima and saddle points are dehned globally just like the lens 
positions, irrespective of which of the contours is the critical curve. 


4. The n-point-mass lens critical curve 
4.1. Critical-curve topology 

The critical curve Zcc is the set of all points in the image plane with zero Jacobian. 
Using the second expression for the Jacobian from equation ffTOj) . we can write the general 
critical-curve condition in the compact form 


9iC = 1 ■ 


( 20 ) 


If we multiply the last expression in equation fITUD by all the denominators, we see that for 
n > 2 the critical curve is a real algebraic curve of order An. The same expression also 
indicates that the zero-Jacobian condition is satisfied if the derivative dzC hes anywhere 
along the unit circle. Hence, the points of the critical curve are the solutions of 


E 






= e-2i<A 


( 21 ) 


where 0 G [0, vr) is a phase parameter fjWittI Il990l ). The normalization of the phase is 
arbitrary; in our choice one gets the single-lens critical curve in the form Zcc = ± \//Ji 


The critical-curve equation fl2Tl) can be multiplied by all the denominators to yield a 
polynomial equation of degree 2n. It follows that for any value of 0 there are up to 2 n 
different solutions of equation (ETj). By varying the phase from 0 to vr the solutions trace 
out 2n continuous critical-curve segments in the image plane. The choice of sign of the 
phase in equation ensures that while tracing any segment in the sense of increasing 0, 
the Jacobian is positive on its right side, negative on its left side. This can be shown by 
differentiating equations (ITOll and (l2lll to obtain 


det J {zq^i d^;) = ±4 dcj ), 


( 22 ) 


where zq is a point on the critical curve, and dz is a tangential shift along the critical curve 
corresponding to phase shift d0. With this notation zo^idz are points offset perpendicularly 
to the right or left, respectively, from the tangential shift dz. 
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The segments connect together in groups of 1 to 2 n to form closed loops, the full 
set of which constitutes the critical curve. This follows from the continuous and periodic 
dependence of equation fl^ on 0 , together with the fact that the critical curve avoids poles 
and complex inhnity, 2 ; 7 ^ Zj,oo. 


From the limitations on the number of segments per loop it follows that the number 
of loops of the critical curve may be as low as one ( as in the ca se of the single lens or the 
intermediate binary lens), but not higher than 2n fjWittIllQQOl) . Individual loops may lie 
outside one another (as in the wide binary lens), or some may he nested in another (as in 
the close binary lens). The number of loops and their mutual position together dehne the 
topology of the critical curve. 


The topology of the critical curve determines the structure of the caustic, which in 
turn determines the number of images as well as the structure of observable light curves in 
microlensing events. The topology depends on the parameters of the n-point-mass lens, i.e., 
on the masses and positions of the individual components. The actual number of relevant 
parameters for n > 1 is reduced from 3 n to 3 n — 4, after accounting for the freedom 
of choice of the coordinate system (i.e., its origin, orientation, and position-measurement 
unit). Mapping the critical-curve topology in this parameter space is a crucial step toward 
understanding the nature of lensing by an u-point-mass lens. 


4.2. Identifying topology changes in parameter space 


In order to map t he critical-curve t o polog y in the parameter s p ace o f the lens, we 
follow the approach of lErdl fc Schneiderl (1199311 and IWitt fc PettersI (Il99311 . and identify 
parameter combinations that correspond to changes in the topology. These combinations 
dehne boundaries of regions with specihc topologies in the parameter space. In order to hnd 
them we hrst note that individual loops of the critical curve are typically disconnected and 
have no self-intersections. Intersections occur only at the boundaries, where the topology 
changes by the splitting or merging of loops. 

The critical curve, viewed as the det J = 0 Jacobian contour, separates regions of 
positive and negative Jacobian in the image plane. Therefore, at a simple intersection point 
on a critical curve where two loops meet, there is a facing pair of positive regions and a 
facing pair of negative regions. At the intersection the Jacobian peaks at zero when crossing 
from one negative region to the other. At the same time it drops there to zero when crossing 
from one posi tive region to the other. Hence, such an intersection must be a saddle point 
of det J (e.g., ISchneider fc Weisslll986l) . In case three loops meet at an intersection point. 
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the Jacobian does change sign bnt nndergoes an inflexion when crossing it. The three-fold 
symmetry forms a monkey saddle on the snrface, as mentioned in § 13.21 Similarly, higher- 
order intersections of the critical cnrve correspond to higher-order saddles on the Jacobian 
snrface. Generally, intersections of any order mnst simnltaneonsly satisfy the saddle-point 
eqnation (IT^ as well as the critical-cnrve eqnation (1^ . 

The transitions in parameter space can be illnstrated on the example of the binary 
lens in FignredJ For lens separations s > 1.861 the critical cnrve consists of two separate 
loops aronnd the lenses. At s ~ 1.861 the critical cnrve passes throngh a saddle point, at 
which the two loops come into contact. For lens separations 1.861 > s > 0.733 the critical 
cnrve consists of a single merged loop snrronnding both lenses. At s ~ 0.733 the critical 
cnrve self-intersects at a symmetric pair of saddle points, splitting off two small loops. For 
lens separations s < 0.733 the critical cnrve consists of an onter loop pins two small loops 
snrronnding the maxima of the Jacobian. Overall, for the given mass ratio the different 
topologies are separated by bonndaries at s ~ 1.861 and s ~ 0.733, at which the critical 
cnrve passes throngh saddle points of the Jacobian. 


As discnssed in § 13.21 the n-point-mass lens may have np to 3 n — 3 different saddle 
points Zsadd- Each of them lies on a Jacobian contonr det J {z) = det J (zsadd)- While in the 
binary lens two of the three saddles lie on the same contonr, for higher n there may be in 
principle np to 3 n — 3 different contonrs dividing the image plane into regions with different 
critical-cnrve topology. Taking into acconnt the scaling properties discnssed in § 13.31 this 
means that any lens conflgnration may nndergo np to 3 n—3 topology changes when changing 
its scale from the “wide” to the “close” limit. By conflgnration we mean here fixed relative 
masses and positions of the components except an absolnte scale. This property can greatly 
simplify the mapping of critical-cnrve t opologies for different n-p oint-mass lens models, as 
demonstrated in the companion paper (jPanek fc Hevrovskvll2015l) . 


5. The n-point-mass lens caustic 

5.1. Caustic structure and cusp number 

The caustic Cc is a characteristic curve of the lens, defined as the set of all points in the 
source plane with infinite point-source amplification. Since the amplification of an image is 
the inverse absolute value of its Jacobian, the caustic must have an image with zero Jacobian. 
Indeed, in addition to regular images the caustic has a degenerate image, the critical curve 

^cc- 
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If we know the critical curve, the caustic can be obtained simply from the lens equation 


Cc ^cc ^ ^ 
i=i 


jjjj 


(23) 


Since the critical curve passes neither through lens positions, nor through complex infinity, 
the simple continuous transformation in equation fl2^ preserves various properties of the 
critical curve. For example, each loop of the critical curve corresponds to a loop of the caustic, 
so that the number of loops is identical. The boundaries in parameter space corresponding to 
changes in critical-curve topology thus also correspond to changes in the number of caustic 
loops. 


Unlike for single or binary lenses, the mapping between the critical curve and the caustic 
given by equation fl25]) for lenses with n > 3 components is not always one-to-one. Different 
points of the critical curve may be mapped onto the same point of the caustic. As a result, 
loops of the caustic may intersect each other, and individual loops may self-intersect. While 
following a path in parameter space, a caustic loop may move from the inside of another 
loop and cross it to the outside without a change in the number of loops. Hence, the mutual 
position of loops does not play the same role as in the critical curve, where the only contact 
between loops occurs at points where the loop number changes. The mutual position of 
caustic loops does play a role in identifying the number of images of the background source, 
which changes when crossing the caustic, as described in § [2l The caustic thus plays the role 
of a boundary between regions of the source plane with different numbers of images. 


Another distinction of the caustic from the critical curve is its lack of smoothness. 
Within our n-point-mass lens model given by equation ([2]), all loops of the caustic have 
cusps, sharp points at which the tangent-vector orientation changes by tt. The smooth parts 
of the caustic between cusps are called folds. The number of cusps of a given caustic loop 
changes only under special local caustic metamorphoses. 


The range of possible metamorphoses expands when increasing the number of lens com¬ 
ponents. The binary lens with n = 2 permits only the simplest beak-to-beak metamorphosis, 
in which two tangent folds reconnect and form two facing cusps. This metamorphosis corre¬ 
sponds to a change in critical-curve topology; the boundaries are thus identical. 

All higher-order metamorphoses occur without any accompanying change in the critical- 
curve topology. For lenses with n = 3 components, two additional metamorphoses may occur. 
In the swallow-tail metamorphosis two cusps and a self-intersection arise from a point along 
a fold. In the butterfly metamorphosis two additional cusps and three self-intersections arise 
from a cusp. 

In the rest of this section we concentrate on identifying the conditions under which 
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caustic metamorphoses occur. In the parameter space they dehne boundaries that form a 
hner subdivision of the critical-curve topology regions discussed in § 14.21 according to the 
number of cusps of the caustic. 


5.2. Cusp condition 


Cusps of the caustic are points at which the caustic tangent changes phase by tt. Since 
the caustic is computed from the critical curve following equation fl2^ . we derive the tangent 
to the caustic from the tangent to the critical curve. 


We start with the gradient of the Jacobian in complex notation, 




d det J . d det J 
dRe[2;] dlm[;2] 


2 dz det J. 


(24) 


The critical curve is a contour of the Jacobian, hence its tangent is perpendicular to the 
gradient. We dehne the tangent to the critical curve as 


7); = i Gz = 2 i dj det J , 


(25) 


evaluated along the critical curve fjWittI Il990l) . In this dehnition the absolute value of 
is equal to the absolute value of the gradient at the given point. This implies that the 
only points that may occur along the critical curve with = 0 are Jacobian saddle points. 
Critical-curve loops meet there for lens parameters at the boundaries studied in § 14.21 


The orientation of is such that the Jacobian increases (i.e., is positive) to its right 
and decreases (i.e., is negative) to its left. For our lens equation ([5]) we may use the complex 
conjugate of the left equality from equation flTT]) to obtain 


T, = -2idzCa|C. 


(26) 


Using the Jacobi matrix J from equation 
from the transformation 

Tc 
Tc 


we compute the tangent to the caustic Ta 


= J 


Tz 

% 


which gives us flWittlll990ll 


T. = T,d,C + f,d,C- 


(27) 

(28) 


Using lens equation (|5]) and equation (|26|) . we obtain 


Tf = 2i [(d-xfglc-d,cdlc] . 


( 29 ) 
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As we trace the critical curve, the right-hand side of equation fl29l) changes continuously, 
since the critical curve avoids the singularities at the lens positions. The only way in which 
the argument of T<^ may change abruptly by tt is when |T^| shrinks to zero. Therefore, cusps 
occur at points along the caustic where = 0. 

Using equation (|29|) we conclude that in the image plane, cusps correspond to points 
along the critical curve which fulhll the cusp condition 

{d,cfdlc = dxdlc. (30) 


5.3. The cusp curve and its properties 

If we multiply equation fl30D by d‘^( {dzCY ^^^1 utilize the critical-curve equation fl20|) . 
we get 

mr = \djc\HdzO^. ( 31 ) 

We note that for Jacobian saddle points d% = 0, so that both sides of equations fl5U]) and 
flST]l are equal to zero. If a saddle lies on the critical curve, it corresponds to cusp images. 
At such points in the image plane loops of the critical curve meet, while in the source plane 
cusps meet in a beak-to-beak metamorphosis, as discussed in § 15.11 

For all other solutions we divide equation flJip by the hrst term on the right-hand side. 
From equation (1211) we know that the last term on the right-hand side lies on the unit circle 
for critical-curve points. We may thus write the result as 

arg = arg (32) 

an equation involving only the arguments (phases) of the powers of lens-equation derivatives. 
An equivalent formulation of equation fl3^ is 

{dlCf=AA{dzCr, (33) 


where A is a real non-negative parameter, as seen from equation fl3T|) . In view of the one 
degree of freedom, equation fl3^ or fl3^ describes a curve in the image plane which intersects 
the critical curve at the positions of cusp images. For brevity we call it the cusp curve. 


If we compute the derivatives of ( using lens equation ([5]), equations fl32|) and fl3^ can 
be converted to 


arg 


-\ 2 


E 

. 1=1 




zjY 


= arg 


E 

. 1=1 


hi 




3 


(34) 
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and 


_^ 

(9- - 


flj 




= A 


n 

z 

Li= 


/ij 


(35) 


Zj)”] [~^ (9 - 2j)' 

The latter equation can be converted to a polynomial of degree 6n — 6, so that varying A 
from 0 to cx) yields up to 6 n — 6 branches of the cusp curve. The cusp curve connects all 
important points of the image plane: the lens positions, Jacobian maxima, Jacobian saddle 
points, and complex inhnity. The value A = 0 corresponds uniquely to the saddle points 
(up to 3 n — 3 different ones), and A = oo uniquely to the maxima (up to 2 n — 2 different 
ones). Complex inhnity corresponds to several roots at A = 1. Similarly, each lens position 
Zj corresponds to several roots at A = l//rj. 


Changing the scale in the image plane {z,Zj) —)■ a{z,Zj) leaves equations fIM)) and fl5S]) 
unchanged, which means the cusp curve is scale-invariant. Therefore, when added to a 
Jacobian-contour plot such as the right panel of Figure [H it will intersect all contours (i.e., 
re-scaled critical curves) at the positions of cusp images. A single plot can thus illustrate the 
number and distribution of cusps not only for the given lens conhguration, but simultaneously 
also for all arbitrarily shrunk or expanded conhgurations. We point out here that the initial 
forms of the cusp condition in equations fl30|l and fl3Tl) are not scale-invariant, hence they 
cannot be used in the same way. 


For illustration we present in the left panel of Figure |2] the cusp curve of the binary lens 
from Figure [T] with pi = 4/5 and Sq = 1. The curve can be constructed from 6 branches 
starting at the saddle points and ending at the maxima, with some branches turning en 
route at the lens positions (two branches through each) or at complex inhnity (two branches 
proceeding to either lens). The intersections of the cusp curve with the critical curve show 
the positions of cusp images, hence the corresponding caustic has six cusps. 


Intersections with the other contours show us the positions of cusp images for re-scaled 
lens conhgurations (in this case for diherent binary separations s). The number of cusps 
changes in this case only through beak-to-beak metamorphoses that occur here for separa¬ 
tions s ~ 1.861 and s ~ 0.733 (see also the discussion in 14.2p . If we increase the separation 
to s > 1.861, two additional cusps appear between the lenses, as shown by the contour 
closest to the lenses. The two loops of the critical curve each have four intersections, hence 
the two caustic loops have four cusps each. On the other hand, if we reduce the separation 
to s < 0.733, two additional pairs of cusps appear outwards from the off-axis saddle points. 
The outer loop of the critical curve has four intersections (its caustic loop has four cusps); 
the two small loops around the maxima have three intersections each (their caustic loops 
have three cusps each). 

The geometry of the cusp curve and the sequence of caustic transitions remains the 
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same for any two-point-mass lens. Adding a third component to the l ens opens np a range 
of diff erent geometries, as illustrated by the examples in § [6] and in iDanek fc Hevrovskv 
(120151) . The cusp curve then becomes a useful guide for understanding the occurrence of 
different caustic metamorphoses, and for keeping track of the cusp number on individual 
caustic loops. 


5.4. Caustic metamorphoses and the morph curve 


When changing the scale of a given lens conhguration, the local appearance of two 
additional cusps other than at beak-to-beak metamorphoses may be detected in the image 
plane by hnding points where the cusp curve is tangent to the local Jacobian contour. For 
a simple tangent point, contours adjacent from one side have no intersection with the cusp 
curve while contours adjacent from the opposite side each have two intersections. Such a 
conhguration corresponds to the swallow-tail metamorphosis. 

If the tangent point occurs at a simple self-intersection point of the cusp curve away 
from saddle points, contours adjacent from one side have one intersection with the cusp 
curve while contours adjacent from the opposite side each have three intersections. Such a 
conhguration corresponds to the butterhy metamorphosis. 


The appearance of additional cusps at a point on the caustic requires stationarity of the 
caustic tangent at the correspondi ng point in the image plane when moving along the 
critical curve le.g.. iPetters et al.ll200l[) . This requirement implies the condition 

{Tz dz + Tz dz) = 0 , (36) 


which must be fulhlled in addition to the critical-curve equation (12Up and the cusp condition 
given by equation (l30l) . 


In the case of our lens equation ([5]) we can substitute from equations (l25l) and (l28l) for 
the tangents Tz and T^, respectively. Equation (15^ can then be written in the somewhat 
ungainly form 

4c9;c[9;ca;c-(4cta;c] =9.ca;c[9.csic-2a.-C(^c9y]. ( 37 ) 

The terms before the brackets indicate that all stationary points of the Jacobian are among 
the solutions of equation (1361) . Here we seek other solutions, such as the swallow-tail and 
butterfly metamorphosis points. In order to hnd them we simplify equation (1371) utilizing 
the critical-curve and cusp conditions. 

We hrst multiply equation (l37p by d^( and substitute for {dzCY the second terms 

on both sides of the equation using the cusp equation (l30l) . Next we multiply the obtained 
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equation by dz(, which cancels out dzC on the left-hand side as well as in both terms on the 
right-hand side, due to the critical-curve equation fl20|) . Now we have 

idler (dio^ - dzC idler die = d-ze mr die - 2 idler idler • ( 38 ) 

The second term on the right-hand side can be added to the hrst on the left-hand side. If 
we divide the equation by | die and transfer the second term on the left-hand side to the 
right-hand side, we are left with the real equation 

3 = Re[2dzeidle)-^dle]. (39) 


Since the imaginary part of the expression in the square brackets is not constrained by the 
equation and thus may be arbitrary, we parameterize it as SiT, where T is a real parameter. 
We may now write 

2a,c53^ = 3(a2o'(i + ir). (40) 

We hrst note that in case the Jacobian has higher-order maxima or higher-order saddles (see 
13.21) . they will solve equation (IT0|) for any T. Maxima never lie on the critical curve, while 
saddles always correspond to topology-changing beak-to-beak metamorphoses (with more 
beaks involved in the case of higher-order saddles). 

Setting these isolated points aside, equation fl59D is a real equation in the complex plane 
depending both on 2 ; and z, while equation fl40|) is a complex parametric equation for ; 2 . Both 
describe the same curve, which we call the morph curve. It passes through the critical-curve 
points corresponding to all caustic-metamorphosis points, such as beak-to-beak, swallow-tail, 
and butterfly. In addition, it passes through the lens positions and complex inhnity. 

If we compute the derivatives of C and C using lens equation ([5]), equations (l3^ and 
(I40|) can be converted to 


1 = Re 


n 




i^j 

iz - Zjf 


n 


E 



(41) 


and 

2 


respectively. Multiplying the latter equation by all denominators converts it to a polynomial 
of degree 6 n — 6, so that varying T from —00 to 00 yields up to 6 n — 6 branches of the 
morph curve. The values T = —00 and T = cx3 correspond uniquely to the saddle points 
(up to 3n — 3 different ones). Each simple saddle point is crossed twice in perpendicular 



E 

.7 = 1 




E 

..7 = 1 




ZjY 


= (i + ir) 


E 

.7 = 1 




ZjY 





















20 


directions by the morph curve, monkey saddles are crossed symmetrically three times, higher- 
order saddles a corresponding higher number times. Turning around any saddle point, the 
arriving morph-curve parameters alternate between T = — cxd and T = cxd. Complex infinity 
and the lens positions each correspond to a pair of roots at T = 0. 

From both versions of the morph-curve equation it is clear that changing the scale in the 
image plane {z,Zj) a{z,Zj) leaves equations flTT|) and fH2|) unchanged, which means that 
even the morph curve is scale-invariant. We can therefore add it to Jacobian-contour plots 
and identify its intersections with the cusp curve away from the lens positions. Intersections 
at the saddle points identify the topology-changing beak-to-beak metamorphoses studied 
earlier. Simple intersections away from saddle points identify swallow-tail metamorphoses. 
Tangent intersections of a single morph-curve branch at a simple cusp-curve self-intersection 
point identify butterfly metamorphoses. For any of these points, the Jacobian contour pass¬ 
ing through it identifies the scale of the lens configuration at which the metamorphosis 
occurs. 


5.5. Distinguishing butterflies from swallow tails 

In order to distinguish between swallow-tail and butterfly metamorphoses, we first check 
the condition for self-intersection points on the cusp curve. An arbitrarily small circle around 
such a point must have more than two intersections with the cusp curve defined by equa¬ 
tion By expanding the equation we obtain the condition 

3 = 2d,C{dlO-^dU. (43) 

Comparing this self-intersection condition with the morph curve equation flT0|) . we find that 
all butterfly points must be morph-curve points with F = 0. From the discussion above 
we know that for F = 0 the morph curve has 2 n roots at the lens positions plus 2 roots at 
complex infinity. It follows that a given lens configuration may have at most 4 n — 8 butterfly 
points. 

Butterfly metamorphoses additionally require the critical curve to be tangent to one of 
the branches of the cusp curve. Expanding the critical-curve equation fl20l) to first order 
around its point zq away from saddles, we get its tangent 

^ ^0 + i (4C <9^0"^ , (44) 

where all derivatives are computed at zq, and the real parameter ei varies along the tangent. 
Similarly, we expand the cusp-curve equation fl32|) to first order around a self-intersection 
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point away from maxima and lens positions. The tangents of its two branches are 



(45) 


where the phase 


^ = argji [3(5y)=> - 5JC(4C)'] [(SzO'SK]"'} . 


(46) 


all derivatives are computed at zq, and the real parameter S 2 varies along the tangent. For 
one of the two perpendicular tangents to be equal to the critical-curve tangent, the hrst- 
order terms of the expansions in equations and fHSjl must have the same argument in 
the complex plane (up to an integer multiple of tt). Setting the arguments equal and using 
equations fl2UD and fl^ we derive the condition 


im|c)y(c)y8;c)-‘] =o. 


(47) 


Any r = 0 morph-curve intersection with the cusp curve has to be checked using this 
tangency condition, to make sure the metamorphosis is a butterfly rather than a pair of 
swallow tails. 

In the rare case when the critical and cusp curves additionally have the same curvature 
at their tangent intersection point, the condition in equation flT7|) is not sufficient. The 
described technique then has to be extended by comparing higher-order expansions of the 
curves. Generally, a butterfly-like transition from one to three intersections occurs if the 
hrst difference between the expansions appears at an even order. Note that the equality 
of each order of the expansions introduces an additional constraint on the lens parameters. 
Such situations would therefore occur in lenses with a sufficiently high dimension of their 
parameter space (i.e., with a higher number of components). 

Finally, we have to check whether we are at a simple self-intersection of the cusp curve, 
since equation flT3|) is valid even for higher-order self-intersections with three or more inter¬ 
secting branches. This hnal condition for the butterfly metamorphosis is 


9.‘03(a.c)-y9.yT. 


(48) 


Violation of this condition would lead to higher-order metamorphoses, which can be analyzed 
in a similar manner. Nevertheless, they do not occur in the triple lens and their exploration 
goes beyond the scope of this paper. 

To summarize, the butterfly metamorphosis occurs when the critical curve intersects the 
cusp curve so that is passes as a tangent through a self-intersection point involving exactly 
two branches of the cusp curve. The conditions on the first through fourth derivatives of 
the lens equation are given by equations fl20|) . fl32l) . flT3l) . flT7|) . and flT8|l . Mapping the 
self-intersection back to the source plane using lens equation ([5]) yields the position of the 
butterfly metamorphosis point on the caustic. 
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6. Lens analysis using cusp and morph curves 


In the right panel of Figure |2] we present the morph curve of the binary lens from 
Figure [T] plotted together with the cusp curve. Their hve mutual intersection points include 
the two lens positions and the three Jacobian saddles. At the lens positions the incoming 
cusp-curve branches alternate with the morph-curve branches. At the (simple) saddle points 
the cusp curve crosses the intersecting morph-curve branches, identifying the occurrence of 
beak-to-beak metamorphoses for correspondingly re-scaled lens conhgurations. The morph 
curve has the same geometry for all two-point-mass lenses, marking the positions of the three 
beak-to-beak metamorphosis-point images along the respective Jacobian contours (re-scaled 
critical curves). 


As mentioned above in § 15.11 swallow-tail and butterfly metamorphoses occur in n- 
point-mass lenses starting from n = 3. We note here that they occur al r eady for a two- 
point-mass lens with non-zero external shear, as shown bv iGrieger et al.l fll989l) . Such a 
model is relevant to the present work as a limiting case of a two-point-mass lens weakly 
perturbed by a third component. In Figure [3] we illustrate the swallow-tail case as seen 
in the image plane of a triple equal-mass lens in a linear conhguration with lens positions 
{zi, 2 : 2 , 2 : 3 } = {—29/30,—2/30, 31/30} s for s = 0.5. Here the scale parameter s is one half 
of the separation of the two outer lenses. 


The Jacobian surface has a generic triple-lens topology with 4 maxima and 6 saddle 
points, which correspond to beak-to-beak metamorphoses as indicated by the intersections 
of the cusp and morph curves in the right panel. In addition, there are 4 simple intersections 
of the cusp and morph curves in the right panel, seen also in the left panel as points where 
the cusp curve is tangent to Jacobian contours (two pairs above and below the axis between 
the maxima). Since they occur away from cusp-curve self-intersections, they correspond to 
swallow-tail metamorphoses. 

We can now directly read off the sequence of gradual caustic changes from Figure [31 as 
we vary the scale of the lens from the wide to the close limit. In the hgure this corresponds 
to progressing from the lowest to the highest Jacobian contour, starting from the small loops 
around the lenses. In the wide limit the caustic has three loops with four cusps each. The 
loops corresponding to the two left bodies merge through a beak-to-beak metamorphosis 
forming a loop with six cusps. Another beak-to-beak metamorphosis merges it with the 
third-body loop, leaving a single 8 -cusped loop. 


The critical curve then reaches the symmetric pair of swallow-tail points closer to the 
lens axis. Two simultaneous swallow-tail metamorphoses add four cusps to form a single 
12-cusped loop. Next comes the symmetric pair of saddle points near the two left maxima. 
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Two simultaneous beak-to-beak metamorphoses split off two 3-cusped loops, leaving a third 
10-cusped caustic loop. The symmetric pair of swallow-tail points further from the lens 
axis changes the third loop by two simultaneous swallow-tail metamorphoses removing four 
cusps, so that it is left with only six cusps. The symmetric pair of saddle points near the two 
right maxima leads to the final pair of simultaneous beak-to-beak metamorphoses, splitting 
off two more 3-cusped loops and thus removing two cusps from the 6-cusped loop. In the 
close limit the caustic has four 3-cusped loops and a single 4-cusped loop. 

In Figure H] we illustrate the butterfly metamorphosis in the image plane of a sim¬ 
ilar triple equal-mass lens, in a linear conhguration but with symmetric lens positions 
{zi,Z 2 ,Z 3 } = { —1,0,1} s for s = 0.5. Even in this case the Jacobian has 4 maxima and 
6 saddle points, which correspond to beak-to-beak metamorphoses indicated by the intersec¬ 
tions of the cusp and morph curves in the right panel. In addition, there are 4 intersections 
along the imaginary axis in the right panel, each with a single branch of the morph curve 
passing through a self-intersection point of the cusp curve. In the left panel we see them as 
points where one branch of the cusp curve at a simple self-intersection point (away from lens 
positions and Jacobian maxima) is tangent to the local Jacobian contour. This confirms the 
intersections correspond to butterfly metamorphoses. 


The sequence of gradual caustic changes is simpler here as we vary the scale of the lens 
from the wide to close limit, due to the higher degree of symmetry of the lens. In the wide 
limit the caustic has three loops around the lens components with four cusps each. All loops 
merge through two simultaneous beak-to-beak metamorphoses forming a loop with eight 
cusps. Two simultaneous butterfly metamorphoses due to the butterfly points closer to the 
real axis add four cusps to form a single 12-cusped loop. Four simultaneous beak-to-beak 
metamorphoses due to the four off-axis saddles then split off four 3-cusped loops, leaving a 
fifth 8-cusped caustic loop. Finally the outer pair of butterfly points causes two simultaneous 
butterfly metamorphoses that remove four cusps from the fifth loop, leading to the close limit 
with four 3-cusped loops and a single 4-cusped loop. 


In Figure [5] we summarize the local behavior of the hrst three elementary caustic meta¬ 


morp hoses occurring in u-point-mass lenses (see e.g., ISchneider et ahl Il992l : iPetters et ah 


20011) . The left, central, and right pairs of columns illustrate the beak-to-beak (vicinity of 


central saddle point from Figure [2]), the swallow-tail (swallow-tail point close to top right 
maximum from Figure [3]), and the butterfly (butterfly point just above central lens from 
Figure Hj) metamorphoses, respectively. The left panel of each pair shows the detail in the 
image plane including Jacobian contours with constant det J spacing. The right panel shows 
the corresponding detail of the caustic in the source plane. Rows correspond to different lens 
scales parameterized by s, as marked in each pair (top: before metamorphosis, middle: at 
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metamorphosis, bottom: after metamorphosis with two new cusps). The image-plane details 
are plotted in units of s, so that the three critical curves in each column correspond to three 
different contours from the same Jacobian contour plot. 


In addition to the three types of local behavior shown in Figure [5l for a triple lens the 
only more complicated metamorphosis occurs at monkey saddles, where three loops of the 
critical curve meet. In the source plane they correspond to beak-to-beak metamorphoses 
with t hree “beaks” meeting at a single point. An example can be found in the companion 
paper (jPanek fc Hevrovskvl 120151) . 


In lenses with n > 3 components higher-order saddles may occur, corresponding to 
a meeting of a higher number of “beaks” in the source plane. In addition, higher-order 
metamorphoses may occur at cusp-curve self-intersections involving three or more cusp- 
curve branches. Their exploration, which may proceed in a similar manner as the butterfly 
analysis in § 15.51 is beyond the scope of this paper. 

As a hnal caveat, note that all three scenarios illustrated in Figure [5] may be complicated 
in special situations. The beak-to-beak scenario may be altered if the cusp curve crosses the 
saddle as a tangent to the critical curve. While this would have no effect on the change in 
critical-curve topology, the caustic may possibly undergo no change in cusp number. In the 
swallow-tail and butterfly scenarios there may be no change in cusp number if the critical and 
cusp curves have the same curvature at the tangent point. The outcome of these situations 
can be studied by expanding both curves around the metamorphosis point and checking 
whether the hrst order at which the curves differ is even or odd. 


All caustic metamorphoses studied here were demonstrated by changing the lens scale, 
which permits the usage of contour plots to visualize the full sequence. However, the meta¬ 
morphoses occur even when changing arbitrary other parameters of the lens. Nevertheless, 
they happen exactly at the points identihed by the presented analysis. All relevant equa¬ 
tions and the sequence of caustic changes remain the same for any crossing of the respective 
parameter-space boundary. 


7. Summary 

In this work we studied n-point-mass lensing by exploring the geometry of the Jacobian, 
viewed as a two-dimensional surface dehned over the image plane. If we plot the Jacobian 
contours in the image plane, the zero-Jacobian contour represents the critical curve of the 
lens. However, as we show in § 13.31 any contour of the Jacobian det J [z) = \ corresponds 
to the critical curve of a lens with component spacing re-scaled by the factor — A. 
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For a lens consisting of a given combination of point masses in a given formation, this 
property permits the simnltaneons study of critical curves for all re-scaled conhgurations, 
from the single-lens close limit to the n-independent-lenses wide limit. A single contour plot 
thus illustrates the lens properties for a full ID cut through the parameter space of the lens, 
instead of just for a single parameter combination. 

We extend the approach further in § [5] by introducing two scale-invariant curves: the 
cusp curve given by equation fl5^ . and the morph curve given by equation (jUj). The cusp 
curve identihes the positions of cusp images along all the Jacobian contours. The morph 
curve pinpoints positions along the cusp curve, at which the number of cusps changes due 
to caustic metamorphoses. 


For illustration, we present several simple applications in § O The overall structure 
of the curves for the binary lens shown in Figure [2] is generic, valid for any mass ratio. 
However, adding a third component already opens up a broader range of possible structures 
and associated lens properties. We use triple-lens examples to illustrate two additional 
caustic metamorphoses, the swallow tail in Figure [3] and the butterfly in Figure 01 Details of 
the image and source planes in Figure [5] demonstrate the characteristic features associated 
with these metamorphoses as well as the simple beak-to-beak metamorphosis. 


In the companion paper flDanek fc Hevrovsk'si l2015l) we use the presented approach 
for a systematic exploration of several simple triple-lens models. However, the concepts 
introduced here can be used also for more general studies. For example, a detailed analysis 
of the structure of critical, cusp, and morph curves would enable a full classihcation of local 
lensing behavior of triple, quadruple, etc. lenses. We conclude that image-plane analysis of 
Jacobian contour plots may prove to be a useful tool for studying and visualizing the nature 
of lensing by n-point-mass lenses. 


We thank the anonymous referee for helpful suggestions. Work on this project was 
supported by Czech Science Foundation grant GACR P209-10-1318 and Charles University 
grant SVV-260089. 
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Fig. 1.— Correspondence of Jacobian contours and critical curves of re-scaled lens con¬ 
figurations, illustrated on a binary lens with fractional masses {/ri,/i 2 } = {4/5,1/5} and 
positions {zi,Z 2 } = (—1/2,1/2}, i.e., separation sq = 1. Left panel: Jacobian surface of 
lens plotted over the image plane z, truncated below at det J = —15, with contours plotted 
for det J e { - 11 . 0 , -2, 0 , 0.5, 0 . 711 , 0.82, 0 . 88 , 0.92, 0 . 945 , 0.962, 0.973}. Right panel; 
same contours with additional innermost det J = —100 plotted in image plane, with the crit¬ 
ical curve plotted in bold. Each contour is identically equal to the critical curve of a binary 
lens with separation and units on axes multiplied by ^1 — det J. Several separation values 
are included as labels for Jacobian contours listed in bold. Crosses: lens components; pluses: 
Jacobian saddle points. The two Jacobian maxima lie within the off-axis critical-curve loop 
sequences. 
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Fig. 2.— Cusp curve (orange) and morph curve (green) of the binary lens from Figure [T] 
The cusp curve intersects Jacobian contours at the positions of cusp images along the critical 
curves of the corresponding re-scaled lens conhgurations. The morph curve intersects the 
cusp curve at the lens positions and at the images of caustic-metamorphosis points for the 
corresponding re-scaled lens conhgurations. For a sequence of contours progressing through 
any such point the number of cusps on the corresponding caustics changes (see also Figure |5]). 
Here the intersections at Jacobian saddle points correspond to beak-to-beak metamorphoses. 
Notation as in Figured! 

















31 



-1 0 1 


Re[z] 



-1 0 1 
Re[z] 


Fig. 3.— Swallow-tail metamorphosis in the image plane: cusp curve (orange) and morph 
curve (green) of an equal-mass triple lens with positions {zi,Z 2 ,Z 3 } = {—29,—2, 31} s/30, 
where s = 0.5. Jacobian contours are plotted outward from the lenses for det J G 
{-1599, -99, -16.6, -5.25, -2.06, -0.69, 0, 0.36, 0.578, 0.698, 0.788, 0.848, 0.891}. The 
cusp-curve branch connecting the two Jacobian maxima above the real axis and its mir¬ 
ror counterpart below the real axis each pass through two points where they are tangent to 
the contours. These points, pinpointed in the right panel as simple intersections with the 
morph curve, correspond to swallow-tail metamorphoses. There are six additional saddle- 
point intersections corresponding to beak-to-beak metamorphoses. Notation as in Figured! 
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Fig. 4.— Butterfly metamorphosis in the image plane: cusp curve (orange) and morph 
curve (green) of an equal-mass triple lens with positions {zi,Z 2 ,Z 3 } = { —1,0,1} s, 
where s = 0.5. Jacobian contours are plotted outward from the lenses for det J G 
{-1599, -99, -16.6, -5.25, -2.06, -0.69, 0, 0.36, 0.578, 0.698, 0.788, 0.848, 0.891}. The 
two cusp-curve branches connecting the Jacobian maxima above the real axis and their 
mirror counterparts below the real axis each pass through a self-intersection point at which 
they are tangent to the contours. These points, found in the right panel as tangent intersec¬ 
tions with the morph curve, correspond to butterfly metamorphoses. There are six additional 
beak-to-beak metamorphoses at saddle-point intersections. Notation as in Figure [H 
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Fig. 5.— First three elementary caustic metamorphoses of the n-point-mass lens. Left 
column: beak-to-beak (detail from Figure [2]); central column: swallow-tail (detail from 
Figure [3]); right column: butterfly metamorphosis (detail from Figure 0]). Left panel in each 
pair: detail of image plane divided by scale parameter s (black: critical curve; gray: other 
Jacobian contours; orange: cusp curve; green: morph curve). Right panel in each pair: 
detail of source plane (red: caustic). Rows correspond to scale parameters s marked in the 
panels. Top: before metamorphosis with n cusps; middle: at metamorphosis; bottom: after 
metamorphosis with n -|- 2 cusps. 














































































